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Abstract 

The breaking rate of an atomic chain stretched at zero temperature by a 
constant force can be calculated in a quasiclassical approximation by finding 
the localized solutions ("bounces") of the equations of classical dynamics in 
imaginary time. We show that this theory is related to the critical cracks 
of stressed solids, because the world lines of the atoms in the chain form a 
two-dimensional crystal, and the bounce is a crack configuration in (unstable) 
mechanical equilibrium. Thus the tunneling time, Action, and breaking rate 
in the limit of small forces are determined by the classical results of Griffith. 
For the limit of large forces we give an exact bounce solution that describes the 
quantum fracture and classical crack close to the limit of mechanical stability. 
This limit can be viewed as a critical phenomenon for which we establish 
a Levanyuk-Ginzburg criterion of weakness of fluctuations, and propose a 
scaling argument for the critical regime. The post-tunneling dynamics is 
understood by the analytic continuation of the bounce solutions to real time. 

PACS numbers: 73.40.G7c, 62.20.MA;, 64.60.Q6 
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I. INTRODUCTION 



A stretched solid breaks instantaneously when the external widening stress exceeds a 
critical value corresponding to the limit of mechanical stability of the system. For smaller 
stresses the fracture is time-delayed and assisted by thermal fluctuations. At sufficiently 
low temperatures there can be quantum effects, and the breaking rate becomes temperature 
independent. This last situation is the main subject of this paper. 

We consider a one-dimensional system (for example, a polymer chain or nanotube) 
stretched by a constant external force applied at the chain ends at zero temperature, where 
quantum tunneling is the only mechanism responsible for the fracture. There have been 
several attempts in the past to study this problem |l| however its many-body nature was 
not taken into account until the pioneering work of Dyakonov |J, who gave qualitatively 
correct results in the limits of weak and large external forces. A more precise analytical ap- 
proach to the problem valid for small external forces has been suggested recently by Levitov 
et al [§. 

We will adopt a general formalism that allows us to look at various limits from a single 
viewpoint. According to Feynman || the tunneling rate per unit length is given by w oc 
| exp( — ^4./^) | 2 where A is the imaginary-time Action calculated along particle trajectories 
connecting the "initial" (stretched) and "final" (broken) states of the chain, h is Planck's 
constant, and the summation is performed over all trajectories. In the quasiclassical limit 
the trajectories contributing most to the sum satisfy the stationarity condition 5A = 0. 
These trajectories describe how the system traverses the classically- forbidden region, going 
from an unbroken chain into a configuration from which even a classical chain released at 
rest would break irreversibly. It is convenient to combine this motion with its time reversal 
to form a localized path called a "bounce" ||. In the quasiclassical approximation the 
tunneling rate is given by 

w = Bexp(-A bounce /h) (1) 

where the amplitude B and action A bounce are calculated for the bounce (there is no factor 
of 2 in the exponential from the squaring because the bounce covers both the past and 
future). We will demonstrate below that the problem of finding bounce solutions is identical 
to the classical problem of finding the equilibrium crack configuration for a two-dimensional 
solid stretched by a constant uniaxial stress - yet another example of the correspondence 
between one-dimensional quantum field theory and two-dimensional statistical mechanics 
@. In view of this equivalence the results of || and || are already implied by the classical 
works of Griffith [§]. 

The present contribution to this field contains two main results. First, in the limit of a 
large applied force close to the limit of mechanical stability we find an exact bounce solution; 
simultaneously this for the first time determines the equilibrium crack configuration in the 
large stress classical limit. Establishing the range of applicability of this solution reveals that 
fluctuations inevitably come into play in the immediate vicinity of the limit of mechanical 
stability. 

Second, we show how the bounce solutions can be used to provide insight into the real 
time post-tunneling dynamics. 



II. GENERAL FORMALISM 



The imaginary-time Action describing a one-dimensional chain stretched by an external 
force p applied to the chain ends has the form 



where p is the linear particle density, x is the spatial coordinate (along the chain length), 
t is the imaginary time coordinate, u(x,t) is the particle displacement field, and dots and 
primes stand for the imaginary time (t), and spatial (x) derivatives. 

The first integral is over all positions and times except for a small region near x — 0, 
where the break will appear. This separates the line into two pieces, so that for all t the field 
u is discontinuous: u(x = +0,t) — u(x = —0,t) = h(t); however, we expect h(t) to be small 
outside a well-defined time interval. We have defined h(t) so that the distance between the 
two parts of the chain is a + h, where a is the equilibrium (p = 0) interparticle spacing in 
the chain. The coupling between the segments of the chain in the tear region is given by the 
potential of cohesive forces V(h) which describes the interaction between the two ends of the 
broken chain for arbitrary separation - that is, it goes beyond the harmonic approximation. 

The properties of the function V(h) can be summarized as follows (Fig. 1). For \h\ <C a 
it obeys Hooke's law: V(h) = pc 2 h 2 /2a, so that the two parts of the chain are joined 
into one elastic medium. For large negative h, the cohesive potential increases without 
bound, V(h — ► — oo) — ► +oo, reflecting the impossibility of indefinite compression. For 
large positive h it approaches a constant, V(h — ► +oo) = 27, which is the work done in 
infinitely separating the halves of the chain; 27 can also be called the bond energy and the 
factor of 2 is introduced to indicate that there is an energy 7 associated with each free end 
of the half- infinite chain pieces. In what follows we will also need to know the properties 
of the cohesive force function G(h) = dV/dh shown schematically on Fig. 2: for small 
h it grows linearly with h,G(h) = pc 2 h/a, reaching a maximal value p c = pc 2 = 7/a at 
h = h c = a, while G(h) decreases to zero as h — > +00. The parameter h c can be called the 
critical bond lengthening, while p c is the limit of mechanical stability of the system - indeed 
only for p < p c does U(h) = V(h) — ph (the total potential energy in the external field) have 
a metastable minimum at h — /ii(Fig.3) corresponding to the equilibrium stretched chain. 

The bounce solution is shown schematically on Fig. 4 using the imaginary time variable 
t as a second (space-like) coordinate. The particle world lines belonging to the two ends of 
the breaking chain deviate significantly from each other for a time interval and then come 
back together. The external force p applied at the chain ends is time- independent, therefore 
it is shown schematically by series of arrows of the same length pointing in the spatial (x) 
direction. 

The Action (0) and the Figures can be interpreted as representing a two-dimensional 
"crystal" of particle world lines subject to a uniaxial stress p, and in this interpretation 
the bounce of Fig. 4 is a critical crack, poised between the small perturbations that can 
spontaneously heal and the large disruptions that lead to fracture. We note in passing that 
our point of view that the crack opening is nonzero everywhere is different from that of many 
classical approaches M , which restrict the break to a finite region completely surrounded by 
an elastic medium. Our treatment contains the standard approach as a special case. 





To see the connection to the classical fracture problem explicitly let us seek extremal 
paths for the Action @ by varying u while keeping the boundary values h(t) fixed. The 
condition 5A/5u = reduces to the Laplace equation 



d 2 u/dt 2 + c 2 u" = (3) 
The solution to (Q) satisfying the boundary conditions u(x = ±0,t) = ±h(t)/2 has the form 

. . p signx f+°° du . . 

u(x,t) = — -xA — / — Muj) exp [iujt - ( \u \\x \/c)\ 

pc z 2 J-oo 2n 

p x r+™ dt'h{t') 

= p^ X + 2^ J-oo it' -t) 2 + x 2 /c 2 { ' 

where the first term describes the stretched chain, the second is due to the break, and the 
Fourier transform h{uj) = J h(t) exp(—iujt)dt has been introduced. Substituting (£|) in @, 
integrating over x, and ignoring overall constants we find the Action in a form that depends 
only on the field h 

A = -^-f h{t)dt\n{\t-t'\/tMt')dt' + [ dt[V(h)-ph) (5) 
2 J\t-t'\>to J 

where D = pc/2ir and to — a /c is a cutoff introduced to prevent singularities for t — > t'. If 
imaginary time is viewed as a space-like variable then the Action (|^) can be interpreted as 
a Hamiltonian describing a "crack" in a two-dimensional "crystal" of particle world lines. 
The crack can be represented by a distribution of fictitious dislocations [[]], and then the 
first term of fl5|) describes the logarithmic interaction of these dislocations with each other 
with an interaction strength D, a "dislocation density" given by — h, and "Burgers vectors" 
pointing along the chain. Similarly the second integral in (JD can be thought of as describing 
the cohesive interaction V(h) between the sides of the "crack" in the presence of the external 
opening field p. A free-energy functional similar to (|^) has been proposed in [llj and (Tl| 
as a field-theoretical starting point of classical fracture mechanics. 

The configuration h(t) for which ([|) is extremal is determined by the condition 5A/5h = 

0: 

D J*~ P^dt? = G(h) - p (6) 

Here and below the singular integral is taken as a principal value. A model of cracks based 
on a singular nonlinear integral equation of the form (§) has been analyzed by Blekherman 
and Indenbom |T2| . 

The sound velocity c sets the upper limit on the value of \h\, therefore the physical 
solutions to @ should satisfy \h(t)/c\ < 1; at the same time the field h(t) can take on 
arbitrarily large values. 

A bounce solution to (||) will have the following form (Figs. 3, 4). For t — » — oo the 
function h(t) starts at the metastable minimum of the total potential U(h),h = hi, then 
makes an excursion past the unstable maximum, h = h 2 , and comes back to h = hi as 
t — > +oo. 



Substituting (]|) into (|j) and subtracting from the result the Action of the stretched 
unbroken chain (h = hi), we find the tunneling Action that goes into the probability ([!]) 

A bounce = f + °° dt[U(h) - U(h) - h - h ^ dU ^) ] (7) 



where we used the result due to Blekherman and Indenbom [TJ] that f_™ dt[G{h) — p] = 



(total zero force along the tear) for the solution to (H) satisfying h(±oo) = hi. 

For large \t\ we can use a harmonic approximation, since the function h(t) is close to its 
asymptotic limit hi. Thus we introduce ip = h — hi and approximate the right-hand side of 
(^) by U"(hi)(fi, valid for \ip\ <C a. Expanding the integrand of ([]) in t'/t and noting that 
h(t) is even in t we find 

However, this is not very interesting, since within the harmonic approximation the tunneling 
Action ([?]) vanishes. 

Further progress can be made by looking separately at the cases of weak and strong 
tearing force where the tunneling Action can be computed in a controlled fashion. 



III. SMALL-FORCE LIMIT, P< P c 

In studying the weak force limit, p <C p c , we can import some ideas from classical fracture 
mechanics. The position of the unstable maximum (h = h 2 ) of the function U(h) (Fig. 3) 
shifts to infinity as p — > 0. However the amplitude of the bounce solution must be larger than 
h 2 , thus implying that over some time the distance between the edges of the tear is much 
bigger than equilibrium interparticle spacing. Since the bounce tail (§) does not contribute 
to the tunneling Action (|?p in quadratic order, we can take h to be nonzero only within a 
time interval — L < t < L(Fig.4). The tunneling time 2L will be determined later. The 
cohesive force, G(h), is negligible for — L < t < L and operative beyond this interval. The 
latter feature implies that there is a strain outside the segment [— L; L] that takes care of 
the assumed h{t) = 0. Then (H) simplifies to 

for — L < t < L. Neglect of the cohesive forces within the tunneling interval implies that 
the edges of the tear are force free, i.e. v!{x = ±0,t) = 0; this condition applied to the 
representation (J|) with the limits of the integration set by — L and L indeed reproduces ([|). 

In the same approximation the total potential U(h) = V(h) — ph can be replaced by 
U (h) = 27 — ph, and (0) simplifies to 

A bounce = f +L 'dt(2 7 -\ph) (10) 



where for p <^ p c we neglected U{hi) + phi/2 = —p 2 a/pc 2 compared to 27. The solution to 
(§) satisfying h(±L) = is the elliptical crack of fracture mechanics [14| : 



Substituting this into fllCf ) and evaluating the integral we find 



A bmmce (L) = A 1 L-p 2 L 2 /AD (12) 

This expression determines the Action parameterized by the time 2L, which is chosen so 
that Af, ounce (L) is extremal: then the Action is stationary with respect to all variations in 
the path connecting broken to unbroken configurations. Eq. (O) has a maximum at 



2L G = l^D/p 2 = 8-fpc/np 2 (13) 

which is identified as the tunneling time. In the context of classical fracture mechanics this 
is called the Griffith criterion || || . Combining ([!]) , (|I2"D , and (fH|) we find that for p<p c 
the fracture rate (per unit time and per unit length of the chain) is given by 

w S (l/cL 2 G )exp{-A bounce (L G )/h} = (p 4 /^D 2 c) exp(-16 7 2 £/^ 2 ) 

= (p 4 /7 2 P 2 c 3 )exp(-87 2 pc/vr^ 2 ) (14) 

where the prefactor B from (|l]) has been estimated by arguing that the break can occur 
anywhere in two-dimensional space-imaginary time: the prefactor \jcL 2 G is the density of 
independent bounces that can fill up two-dimensional space- imaginary time (to be indepen- 
dent, two bounces must be separated in time by more than L G and in space by more than 
cL G ). 

The effective mass involved in the tunneling can be estimated from an argument parallel 
to that of Dyakonov . The representation ([|) implies that the effect of an inhomogeneity 
of size L G in the t-direction perturbs the picture of the world lines over a distance of order 
cL G (vn the x- direction) containing about cL G /a atoms. Each of them has a mass m = pa, 
therefore the effective mass M involved in the tunneling is found to be 

M ^ m(cL G /a) ^ m(pc 2 -f/ap 2 ) = m(p c /p) 2 (15) 



The tunneling time 2L G flTB|), the maximal width of the tear [see Eq(O)], h(0) = 
{p/ttD)L g = 87/717), and the tunneling mass (|i~5"|) all diverge as p — > implying that fracture 
through tunneling has many-body nature. 

Here as in fracture mechanics this approach breaks down in the vicinity of t = ±L where 
the derivative h(t) diverges. There is also an associated stress singularity (which we will not 
discuss here) near t = ±L. These defects can be remedied by taking proper account of the 
cohesive forces jl4| with the general conclusion that all the scaling dependencies predicted 



in the framework of the Griffith approximation are correct and that the numerical factors 
entering Eqs.([L2|)- (|14D are accurate in the limit p p c \ these comments apply equally well 
to the quantum fracture problem. 

For example, the range of applicability of Eq([TT|) can be found by requiring that both 
\h\ < c and h ^> a with the conclusion that independently of the value of the tearing force 
p, (|ll|) can be trusted outside a small region - a dozen microscopic time spans of order a/c 
in the vicinity of t = ±L G . 



For p <C p c the tail of the true solution (§) becomes 



<p(t - ±00) = ^0 = a(p/p c )(L G A) 2 (16) 

which is generally small since p <C p c . 

The analysis given above shows that for weak breaking forces (p <C p c ) the true bounce 
solution is approximated quite well by EqflTTD for \t\ < L, and by Eq. (|16"D for \t\ > L, 
outside of the vicinity of \t\ = L = Lq- The determination of the size of the transient region 
(expected to be a dozen microscopic time scales a/c) and the precise functional form h(t) 
inside it require a detailed knowledge of the cohesive potential V(h); it should not affect the 
main conclusions (P~5|) - (P~5j) . 

Apart from numerical and preexponential factors, the results ([13]) and ([14]) were first 



found by Dyakonov j| using heuristic arguments based on the tunneling properties of a 
particle of effective mass ([15]) passing through a triangular potential barrier. A conformal 
mapping technique similar to what is employed in fracture mechanics was used by 
Levitov et al ||; this reproduced Eq([Ll]) and (0), the numerical factors entering in Eq.(p~3"P 
and the exponential of Eq.(|l4]). However, the link to classical fracture mechanics has not 
been previously noticed. 

Experimental verification of the weak tension results will be hindered by the exponen- 
tially large lifetime of the chain [see Eq.fJTJJ)]. Therefore the case of large tearing force (p close 
to p c ), where the lifetime is shorter, seems more important from the practical standpoint. 



IV. VICINITY OF THE LIMIT OF MECHANICAL STABILITY, AP < P c 

For Ap = p c — p <C p c the tunneling still has a collective nature, because the limit of 



mechanical stability, p = p c , is formally similar to the spinodal of magnetic systems [10 
pTfl . in this limit the tunneling involves the vicinity of h — /i c (Fig.2) and the bounce 
amplitude is very small. Regardless of the underlying interactions the cohesive force can be 
approximated here as G{h) = p c — b(h — h c ) 2 where b = —^d 3 V(h = h c )/dh 3 . The zeros of 
G(h) — p = Ap — b(h — h c ) 2 determine the positions of the metastable minimum (hi) and 
unstable maximum (^2) of the total potential energy function U(h) (see Fig. 3): 



h h2 = h c T\/Ap/b (17) 
Expressing G(h) — p in terms of if = h — h\ we find 

G(h) -p = 2sjkpb<p - b V 2 (18) 
Similarly the potential energy function U(h) = V(h) — ph can be written as 

U(h) = Uih) + ^Apb v 2 - (19) 
and Eqs. (|]) and ( |T8"P give the equation for the bounce: 



-DO 



0(f) dt' 



D I ^y-r = y A P bi P - V (20) 

-00 v — X 



Close to the limit of mechanical stability the properties of strongly stretched unbroken 
chain substantially deviate from those of the equilibrium one: the linear mass density p is 
significantly reduced, and the lattice looses its stiffness at p — p c . The latter means that 
both the sound velocity c, and the parameter D = pc/2n entering fl20|) vanish as Ap — > 0. 
The functional dependencies can be recovered by comparing the harmonic terms of ([19]) and 
(@), y/Apb = pc 2 /a, which gives us 

D = (m 2 bAp) l '\ c £ aibAp/m 2 ) 1 ^ (21) 

The solution to Eq.(^) is 

<P® = ( 22 ) 

where 



$ = AJAp/b (23) 



£ = - {mVbAp)Ui m 

Combining fll9|) and (|22|) -(p4l) we see that the general "sum rule" (§) is indeed satisfied. 
The value of the bounce amplitude ip(t = 0) = $ is small close to p = p c ; at the same time 



h(0) = $ + = h c + 3y Ap/b = h 2 + 2^J Ap/b > h 2 , i.e. as expected for the bounce solution, 
the function h(t) goes past h 2 for some time. 

The divergent time scale £ (|24]) can be interpreted as a tunneling time, and the effective 
tunneling mass M can be estimated by replacing Lq by £ in (|15|): M = m(c£/a) = m. We 
see that as Ap —>■ the effective mass involved in tunneling is comparable to the single 
particle mass in agreement with Dyakonov 0]. 

Substituting the potential energy function ([Tj]) into (^) with the bounce solution given 
by (E2|)-(P|), and computing the integral we find for the tunneling Action 



A bounce = 2n 2 DAp/b = (m 2 /& 3 ) 1/4 (Ap) 5 / 4 (25) 

which in view of ([!]) implies that for Ap -C p c the tunneling rate is given by 

w ^ (l/c£ 2 ) exp(-27i 2 DAp/bh) = (6Ap/mV) 1/4 exp[-const{m 2 /tf) l/ \Apf /A /%] (26) 

where the prefactor is taken from (0), Lq is replaced by £ [EqQ2~3])], and the dimensionless 
constant under the exponential is of order unity. 

Our results (|2~5"D and (|26"D in the limit Ap -C p c reproduce those of Dyakonov [|J even 
though he could not find the exact bounce solution (|22|)-(|2~i|). Even in the absence of the 



exact solution it can be seen that Eq(pU|) has a solution of the form cp(t) = 4y Ap/fef(i/£) 
where £ is given by (|4]) and f(z) is a scaling function obeying an integral equation that has 
no free parameters at all. This identifies £ oc (Ap) -1 / 4 with the tunneling time, and then the 
tunneling Action is estimated as proportional to (Ap) 3 / 2 (Ap) -1 / 4 = (Ap) 5 / 4 in agreement 
with the calculation that led to (1251). 



In view of the equivalence between the bounce solutions and equilibrium cracks of the 
corresponding two-dimensional classical problem, Eqs. (p2|)-(24) also give the equilibrium 
crack configuration in the limit of large stresses close to the limit of mechanical stability p c . 
In this correspondence the imaginary time t becomes a spatial coordinate along the crack, 
D is an elastic constant |J [still vanishing according to (^)], £ is the effective length of 
the crack, and the parameter m will have a meaning of a typical binding energy divided 
by interparticle spacing. The rate of the crack nucleation will be given by fl2"E| ) with the 
physical temperature T substituting for the Planck's constant h. 

The classical problem of finding the equilibrium crack configuration in the large stress 



limit close to the stability threshold was previously considered in | TT| . The conclusions 
of that work differ from ours for two reasons. First, it was assumed that the long-range 
nature of the interaction between the different parts of the crack profile through the bulk 
of crystal can be replaced by a finite-range interaction; as a result the equation for the 
profile is differential rather than integral [see (EO)]. However, this is certainly not true for 
materials obeying conventional elasticity theory. Second, it was assumed that for small 
(p, G{h) — p = Ap — bip 2 , instead of our EqflTSp. But this cannot be correct, because for 
if = 0(metastable equilibrium) the cohesive force G(h) balances the external force p and 
therefore G(h) — p must vanish. 



V. LEVANYUK-GINZBURG CRITERION AND CRITICAL REGION 



According to ( [25]) the tunneling Action vanishes upon approach to the limit of mechanical 
stability p = p c . However, the results of the previous Sections were based on a quasiclassical 
picture which is only valid for Af, ounce /h ^> 1, i.e. when the exponential factor of (|]) 
dominates the breaking rate. This condition imposes a range of validity for our results 

Ap/Pc > (b 3 h 4 /m 2 p 5 c f 5 (27) 

that eliminates from consideration the immediate vicinity of the threshold of mechanical 



stability p = p c . This condition is analogous to the Levanyuk-Ginzburg criterion [15| of 
weakness for thermal fluctuations in the theory of critical phenomena. 

The results (PH)-(p6|) were also obtained under the assumption that Ap/p c <C lfonly then 
is the expansion in powers of <p used in ( fL8|) - (|20|) accurate]. This will be consistent with 
(0) only if 

Gi = (b 3 h A /m 2 p 5 c )^ 5 < 1 (28) 

The dimensionless quantity Gi is a property of a given substance and is analogous to the 
Ginzburg parameter of the theory of critical phenomena. In the present context the require- 
ment Gi 1 can be viewed as the condition of weakness of quantum fluctuations. 

The inequalities (p^) and ( p8"D set the range of applicability of the theory of Section 
IV. The immediate vicinity of p = p c corresponding to the reversed sign in the inequality 
(fZ7|) can be called the fluctuational region; quantum fluctuations play here a dominant role 
and the quasiclassical approximation is insufficient. In the case of the reversed sign in the 
inequality (EBT), Gi 1, the theory of Section IV has no range of applicability at all. 



For the case that one of the conditions (p?D, ( p28|) is broken, our analytic approach fails, 
but some understanding can be gained with the aid of a scaling theory constructed in 
analogy with the theory of critical phenomena |nj. Assume that in the critical regime there 



is a single independent divergent time scale £, the tunneling time, and all other physical 
quantities can be expressed in terms of £. According to flT9|), the barrier for tunneling can 
be estimated as (Ap) 3 / 2 /^ 1 / 2 , and therefore the tunneling Action behaves as £(Ap) 3 / 2 /6 1//2 . 
In the critical region the tunneling Action should be of the order of Planck's constant H, so 
that 

f = hb 1/2 /(Apf 2 (29) 

We also expect that the tunneling mass is still of order m which implies c = a/£ for the 
critical behavior of the sound velocity. These conjectures are made plausible by two mild 
self-consistency checks: (i) the divergence in ( p9[ ) is stronger than (|24|)) as expected for the 



critical regime; (ii) the dependencies ( p9|) and fl24]) match each other on the border of the 
critical regime defined by (p7p. 

The conditions (|27|) and (p8|) have the classical analogs that the results of Section IV are 
applicable to the classical fracture of solids whenever the following inequalities are satisfied 

Ap/ Pc > (6 3 T 4 /m 2 ^) 1/s (30) 

Gi = (b 3 T A /m 2 p 5 c ) 1/5 < 1 (31) 
Now the constraints (|30|) and (|3l|) describe the conditions for thermal fluctuations to be 



sufficiently weak so that the critical crack described by is relevant for the case of 

large applied stresses. 

Inside the critical region, Ap/p c <C (6 3 T 4 /m 2 ^) 1 / 5 , thermal fluctuations play a dominant 
role and the mean-field-like treatment of Section IV is invalid. If the Ginzburg parameter 
(|3l|) is not small compared to unity, the theory of Section IV does not apply. 

In the critical regime the analog of (|29|), 

f £ Tb 1/2 /(Ap) 3/2 , (32) 

gives the dependence of the critical crack length £ on the distance to the instability threshold 
Ap. 



VI. REAL TIME DYNAMICS 

Coleman § notes that the analytic continuation of the bounce to real time describes the 
evolution of the system after the tunneling takes place. Introducing the physical time r = 
-it into (|TTD we find that for a small tearing force, p <C p c , the separation between the pieces 
of the chain evolves with time as 

h(r) = 2 ^L g + t* (33) 
pc v 

where L G is the true tunneling time (|13"D. Equation (|33D meets our physical expectations: 
after penetrating through the barrier at t = it = 0, the end of the chain starts moving 



classically from rest from the escape position h(0) = 2pLc/pc = 87/717} which is the turning 
point of the classical motion. The result ( |33|) is also precise analytically as in getting fll!]) 
the potential energy function U(h) was approximated by 2'j—ph which is even more accurate 
in the classically allowed region than under the barrier. 

For small times (t < L), E'g.(|3"3"D predicts a motion with a constant acceleration 2p/pcLc 
that can be understood as follows. Immediately after the break each point of the chain 
excepting a finite segment near the edge will be in mechanical equilibrium. Therefore the 
whole unbalanced force p is applied to the edge segment. The mass of this region is the 
same as the tunneling mass (|HjJ), M = pcLc, and its acceleration due to the external force 
p is of order p/M = p/pcLc in agreement with fl3"3"|). This motion breaks the force balance 
of the region next to the moving edge segment causing the propagation of acceleration away 
from the tear. 

To understand the dynamics inside the chain we use the expression |14| Q for the dis- 
placement field u(x,t) corresponding to the solution (|TT|): 

u ( x , t) = -^ReUx + ictf + (cL G ) 2 ] 1/2 (34) 
pc 

Analytically continuing this to the real time r = —it, we find 



u(x, r) = ^— 2 signx ([(x - cr) 2 + {cL G ) 2 } 1 ' 2 + [(x + cr) 2 + {cL G ) 2 } 1 ' 2 ) (35) 

This expression is a solution to the wave equation d 2 u/dr 2 — c 2 d 2 u/dx 2 = in the two 
parts of the broken chain satisfying the conditions u(±0, r) = ±h(r)/2 and u'(±0,r) = 0. 
It describes sound waves propagating away from the break and leaving behind virtually 
unstrained material moving at constant speed. Not much is happening for r < Lq\ however, 
for r ^> Lq the following picture emerges (by symmetry we only need to consider the x > 
segment). For x cr the strain u' is small and u(x, r) = pr/pc, whereas for x 3> cr, du/dr 
is small and u(x, r) = pxj pc 2 : a moving unstrained region encounters a stretched immobile 
region at x — cr. In a transition region having the width of order cLq the matter is beginning 
to move to become part of the moving region. The unstrained region set into motion has 
mass that grows linearly with time [M(r) = mcr/a = per] and moves at the constant speed 
v = p/pc; the momentum M{r)v = pr is growing linearly in time in response to the constant 
applied force p. Figure 5 shows how this soliton-like wave traveling away from the break is 
presented in terms of particle trajectories. The break occurred at x = 0, r = 0, and the 
Figure was constructed by adding u(x,t) [given by Eq. (|35"D l to x for each r. 

In the large force limit, Ap <C p c , we analytically continue the bounce solution (|2T|) to 
find that the distance between the tearing edges of the chain evolves with time as 



h{r) =h l + $7^^ (36) 



Similar to the case of small tearing force, p p c , the initial stage of the motion, r 
described by Eq(p6|) can be interpreted as a constant acceleration of the tunneling mass m 
by the unbalanced force Ap. However the range of applicability of (^) is rather narrow - 
Eqs(|22|)-(|24|) [and thus (|3~6"D1 were obtained assuming the expansions (0)) and ([19]) to be 
valid. At the same time the distance between the tearing edges grows with time according to 



(|36|), thus implying that as time progresses, h{r)—h\ becomes comparable to the equilibrium 
interparticle spacing a and the expansions (|18D and (|19D cease to be accurate. Therefore (36) 



is valid only for r <C £. In the same approximation the motion in the chain can be obtained 
by substituting (^) in (f|), evaluating the integral, and doing analytic continuation to the 
real time r = —it: 

p c signx ( M/c + £ \ , , 

u {x,t ) = — x + — — [hi + <5>£ ,7 —pr 2 2 37 

pc 2 2 \ [|x|/c + 4j l -r l j 

where the sound velocity c is given by (0). Beyond the range r £, the dynamics is 
expected to be described by a theory qualitatively similar to the one given above for the 
case p -C p c with p c and £ substituting for p and L G . 
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FIGURES 
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FIG. 2. The cohesive force G(h) = dV/dh. It reaches a maximal value p c at the critical bond 
lengthening h c which corresponds to the limit of the mechanical stability of the system. 
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FIG. 3. The total potential energy U(h) = V(h) — ph. In the presence of an external tearing 
force p < p c the total potential energy has a metastable minimum at hi describing the untorn 
stretched chain and an unstable maximum at hi , which is the position of the peak of the tunneling 
barrier. 




FIG. 4. The two-dimensional "crystal" of particle world lines torn by a constant external force 
p traversing the imaginary time (t) direction. The tear takes place at x = 0, the distance between 
the torn edges is h + a, the tunneling time is 2L, and the system penetrates through the barrier at 
t = 0. 
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FIG. 5. The collection of particle trajectories constituting a breaking one-dimensional solid 
torn by a constant external force p. 



